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ABSTRACT 

If the first (PopIII) stars were very massive, their final fate is to collapse into 
very massive black holes. Once a proto-black hole has formed into the stellar 
core, accretion continues through a disk. It is widely accepted, although not con- 
firmed, that magnetic fields drive an energetic jet which produces a burst of TeV 
neutrinos by photon-meson interaction, and eventually breaks out of the stellar 
envelope appearing as a Gamma Ray Burst (GRB). Based on recent numerical 
simulations and neutrino emission models, we predict the expected neutrino dif- 
fuse flux from these PopIII GRBs and compare it with the capabilities of present 
and planned detectors as AMANDA and IceCube. If beamed into 1% of the sky, 
we find that the rate of PopIII GRBs is < 4 x 10^ yr^^. High energy neutrinos 
from PopIII GRBs could dominate the overall flux in two energy bands [10^ — 10^] 
GeV and [10^ — 10^] GeV of neutrino telescopes. The enhanced sensitivities of 
forthcoming detectors in the high-energy band (AMANDA-II, IceCube) will pro- 
vide a fundamental insight on the characteristic explosion energies of PopIII 
GRBs and will constitute a unique probe of the the Initial Mass Function (IMF) 
of the first stars and of the redshift zj marking the metallicity-driven transition 
from a top-heavy to a normal IMF. The current upper limit set by AMANDA- 
BIO implies that such transition must have occurred not later than Zf = 9.2 for 
the most plausible jet energies. Based on such results, we speculate that PopIII 
GRBs, if not chocked, could be associated with a new class of events detected 
by BeppoSax, the Fast X-ray Transient (FXTs), which are bright X-ray sources, 
with peak energies in the 2 — 10 keV band and durations between 10 — 200 s. 

Subject headings: stars: early-type - gamma rays: bursts- neutrinos - black holes 
- cosmology: theory 

1. Introduction 

One of the most challenging problems in modern cosmology is the understanding of the 
first episodes of star formation in the Universe. A growing body of theoretical work has been 
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devoted in the past few years to this subject (Rees 1976; Rees & Ostriker 1977; Silk 1977, 
1983; Haiman, Thoiil & Loeb 1996; Uchara et al. 1996). These investigations are aimed 
at the identification of the characteristic mass scale of the first stars (usually referred to as 
PopllI stars), how this relates to the physical conditions of the gas in the first collapsed 
objects and how it affects subsequent structure formation in the Universe. 

Numerical simulations (Abel et al. 1998; Nakamura & Umemura 1999; Bromm, Coppi & 
Larson 1999, 2001; Abel, Bryan & Norman 2000; Bromm et al. 2001; Ripamonti et al 2001) 
based on hierarchical scenarios of structure formation and/or detailed stellar collapse models 
have shown that the typical mass scale for the first collapsed clumps of primordial gas is 
pa IO^Mq, which corresponds to the Jeans mass set by molecular hydrogen cooling. 

The ultimate nature of the stars that form out of these clumps critically depends on the 
physical conditions of the gas as the evolution pushes density to higher values. Preliminary 
studies (Omukai & Nishi 1998; Nakamura & Umemura 1999; Bromm et al. 2001; Ripamonti 
et al. 2001) show that as long as the metallicity is below some critical value (typically 

~ IO^^Zq) the first clumps have httle tendency to further fragment, as also expected 
from a number of physical arguments (see e.g. Schneider et al. 2001). 

Tentative evidences for an early top-heavy initial mass function (IMF) are provided by a 
number of observations (see Hernandez & Ferrara 2000 and references therein). For instance, 
a comparison between the observed number of metal poor stars with the one predicted 
by cosmological models, implies that the stellar characteristic mass sharply increases with 
redshift. Furthermore, the intracluster medium (ICM) metal abundances measured from 
Chandra and XMM spectral data are higher than expected from the enrichment by standard 
IMF supernova (SN) yields in cluster galaxy members, indicative of a top-heavy early IMF. 
Finally, the observed abundance anomalies {e.g. oxygen) in the ICM can be explained by an 
early generation of PopIII SNe (Loewenstein 2001). 

These issues, highly suggestive of a top-heavy early star formation, have recently moti- 
vated a series of numerical investigations of the nucleosynthesis and final fate of metal-free 
massive stars (Heger & Woosley 2001; Fryer, Woosley & Heger 2001). Stars with masses 
in the range MOM© < M < 26OM0 undergo electron-positron pair instability (SN^-y) and 
end up in a giant, nuclear-powered explosion, leaving no remnant and enriching the am- 
bient medium with their nucleosynthetic products. The kinetic energy released during the 
thermonuclear explosions powered by pair instability are fa 10^ times larger than those of 
ordinary Type II SNe. This might cause the interaction with the circumstellar medium to 
be as strong as predicted for hypernovac (Woosley & Weaver 1982). These explosions do not 
lead to the ejection of strongly relativistic matter and therefore cannot power a Gamma-Ray 
Burst (GRB, Fryer, Woosley & Heger 2001). 
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However, if stars have masses M > 260Mq, photodisintegration instability is encoun- 
tered before explosive nuclear burning can reverse the implosion and the stars collapse to 
Very Massive Black Holes (VMBHs), swallowing virtually all previously produced heavy ele- 
ments. These stars are likely to be rapidly rotating and the estimated angular momentum is 
sufficient to delay black hole formation (Fryer, Woosley & Heger 2001). Once a proto-black 
hole has formed into the core, accretion continues through a disk at a rate which can be as 
large as 1 — 1OM0S~^. It is widely accepted, although not confirmed, that magnetic fields 
might drive an energetic jet which can produce a strong GRB through the interaction with 
surrounding gas. 

In this scenario, the energetic jets generated by GRBs engines produce, by photon-meson 
interaction, a burst of TeV neutrinos while propagating in the stellar envelope (Meszaros & 
Waxman 2001). It is widely assumed that if the progenitors of GRBs are massive, collapsing 
stars (the so called "collapsar", Woosley 1993, Paczyhski 1998) the shocks producing the 
7-rays occur after the relativistic jet has emerged from the stellar envelope. For a significant 
fraction of collapsars, the jet may be unable to punch through the stellar envelope (Mac- 
Fadyen, Woosley & Heger 2001). However, the TeV neutrino signal from such "chocked" 
jets should be similar to that from jets which do break through the stellar envelope, leading 
to observable GRBs. Detecting the neutrino fiux from individual collapsars may provide a 
direct evidence of such 7-ray-dark collapses (Meszaros & Waxman 2001). 

In this paper, we compute the diffuse flux of high energy neutrinos emitted by PopIII 
GRBs at high redshifts (5.4 < z < 30). Hereafter, by PopIII GRBs we indicate PopIII stars 
with masses > 26OM0 which collapse to a VMBHs after a transient phase of mass accretion. 

Following Schneider et al. 2001, we assume that a fraction (1 — /^^) of the first stars 
collapse to VMBHs powering a relativistic jet that generates a fiux of high energy neutrinos. 
Using the neutrino emission model proposed by Meszaros & Waxman (2001) and integrating 
over the source rate throughout the universe, we compute the diffuse neutrino flux from 
PopIII GRBs. 

The sensitivities of present (AMANDAQ, Andres et al. 2001) and forthcoming neutrino 
telescopes (such as AMANDA-II, Barwick 2001, IceCubeQ, Spiering 2001 and other km-scale 
detectors, Halzen 2001) enable to derive important constraints on the flrst episodes of star 
formation, such as the relative number of SN^^ and VMBHs {i.e. PopIII GRBs), metal 
enrichment and the nature of dark matter. 



^ littp : // ajnanda . berkeley . edu/ ainanda/ amanda . html 
^ tittp : / / www . ssec . wise . edu/a3ri/ icecube/ 
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The paper is organized as follows: in Section 2 we describe the proposed scenario for 
the formation of PopIII stars and we compute the expected rate of GRBs. In Section 3 
we determine the energy spectrum of TeV neutrinos emitted in individual Poplll GRBs. 
In Section 4 we compute the expected diffuse flux from PopIII GRBs. In Section 5 we 
explore the parameter space making detailed predictions for the detectability with present 
and forthcoming neutrino telescopes. In Section 6 we illustrate the implications of the 
proposed scenario for PopIII GRBs. Finally, in Section 7 we summarize our main results. 



2. Population III stars as GRB progenitors 

In this Section we discuss the formation and evolution of PopIII stars. We work 
within the paradigm of hierarchical cold dark matter (CDM) models for structure formation, 
wherein dark matter halos collapse and the baryons in them condense, cool and eventually 
form stars. We adopt a cluster-normalized ACDM cosmological model with the following 
parameters: Qm — 0.3, Qa — 0.7, h — 0.65 and Qs^^ — 0.019. 

If a fraction of the first stars arc progenitors of VMBHs and can power a GRBs, we 
expect a large number of TeV neutrinos emission episodes to have occurred throughout 
the Universe. In order to estimate the diffuse flux contributed by this large ensemble of 
uncorrelated neutrino bursts, we need to compute the expected number of sources as a 
function of redshift. Following Schneider et al. 2001, to each halo of total mass M, we 
associate a baryonic mass given by {D,B/fiM)M. Here we consider only halos within which 
gas can cool and form stars. In agreement with numerical simulations (Bromm et al. 2001), 
we assume that only a fraction j\ = 1/2 of the baryonic gas turn into stars, the rest remaining 
in diffuse form. The relative fraction of SN-y-y and VMBH progenitors are parametrized as 
follows: 

M^T = At V = /tt^*, (1) 

where M^^ (M,) is the total object mass which ends up in SN^^ (VMBHs) and is the 
mass processed into stars. Thus, only a fraction (1 — /^^) of the formed stars might generate 
a relativistic jet during the VMBH collapse, leading to a TeV neutrino burst. The remaining 
fraction f^^ contribute to the metal enrichment of the intergalactic medium (IGM) and is 
responsible for the transition to a standard fragmentation mode, hence IMF, which occurs 
when the average metallicity is above the critical value < Zcr > ~ 10~'^Zq at redshift Zf. 
Using the SN^^ metal yields given by Heger & Woosley (2001) (see also Umeda & Nomoto 
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2002) and assuming that metals pollute uniformly the IGM, Schneider et al. (2001) have 
inferred the transition redshift, 5.4 < zj < 11, from a top-heavy to a standard IMF as a 
function of /^^. These parameters directly affect the number of PopIII GRBs (or VMBHs) 
predicted by the model, as they define the duration and efficiency of massive PopIII star 
formation epoch. Here we are interested to redshifts z ^ zj within a range 5.4 < 2; < 30. If 
^ 1, i.e. all PopIII stars are SN^^, metal enrichment is very efficient and transition to 
a normal IMF occurs at Zf ^ 11. Conversely, as the relative fraction of SN^-y to VMBHs 
decreases, massive PopIII star formation can proceed down to a minimum redshift Zf = 5.4 
which corresponds to = 8 x 10^'^. At this epoch, the amount of VMBHs formed 

can account for the entire baryonic dark matter content of galaxy halos (see Schneider et 
al. 2001). 

As our reference model PopIII GRB progenitor, we consider a SOOM© star. This choice is 
motivated by the availability of detailed and complete evolutionary models of zero-metallicity 
stars of this mass (Fryer, Woosley & Heger 2001). These simulations include angular mo- 
mentum transport through all burning stages until complete collapse to VMBHs occurs. 

We compute the star formation rate of PopIII stars as, 

m[z) = —l dMn{M,z)M^ (2) 

at JMn,U^) 

where n{M, z) is the comoving number density of halos per unit mass predicted by the Press- 
Schechter formalism, M^, is given by eq. |l| and the integration is performed from M^[^{z) 
which is the minimum mass that can cool within a Hubble time at the specified formation 
redshift z, i.e. tcooi(^5 z) ^ tn (see Ciardi et al. 2000). Therefore, we assume that all gas in 
halos with M > M^ia{z) is able to form stars with an efficiency f\. 

The rate of PopIII GRBs produced by SOOMq progenitor stars is 



R(z) = ^ ^ / dz-^^^, (3) 

^ ' 3OOM0 Jzf {l + z) dz' ^ ' 

where dV is the comoving volume element, 

dV = Atc (^-|-) r{zY T{Qm, ^a, z) dz, (4) 
Hq = 100 /ikms~^Mpc~^ is the Hubble constant, r(z) is the comoving distance and 

J^{flM, ^A, Z) = ^= (5) 

'{i + zy{i + nMz)-z{2 + z) 



The results are shown in Fig. 1 (solid lines). The maximum and minimum rates correspond 
to the extremes values for /^^ (or equivalently, zj) which define the range of signals that 
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might be detected by present and forthcoming neutrino telescopes: Z™"^ and f-y-^ = 0.98 
{zf = 10). As seen, the total rates range from 4 x 10^ to 8 x 10^ events/yr. If the jets are 
beamed into 1 % of the sky, the effective rate will range between 4 x 10^ to 8 x 10^ events/yr. 

The dotted lines in Fig. 1 represent the predicted duty cycle of the signal. This is defined 
as the ratio of the typical duration of a single burst, tj ~ 10(1 + z) s (see next Section), to 
the mean time interval between two successive bursts [1/-R(z)]. 

The total duty cycle is found to be > 1 but, if zj ^ 10, the effective duty cycle might 
be significantly reduced by beaming. This means that the diffuse flux is stationary but 
might not be continuous in time, being characterized by a sequence of individual bursts 
whose average separation depends on the rate. Integrating over a typical observation time 
of ~ 1 yr, the expected number of bursts is > 10^ depending on the assumed values for 
/^^ and on the beaming factor (see Fig. 0). Therefore, we assume that the average diffuse 
flux provides a good representation of the signal averaged over an observation time of 1 yr^ 
Signal detection from individual bursts by forthcoming km-scale neutrino telescopes will be 
discussed in Section 6. 



3. The Emission Model 

In this Section we derive the expected neutrino flux contributed by a single source, i.e. 
a zero-metallicity, rotating, 300Mq star which collapses to a VMBH. 

The simulations carried out by Fryer, Woosley & Heger (2001) show that the structure 
of the pre-coUapse star is that of a I8OM0 He-core with radius 6 x 10^^ cm surrounded 
by a large H-envelope which extends up to ~ 2 x 10^*^ cm. After helium burning, the 
core encounters electron-positron pair instability and collapses igniting explosive oxygen and 
silicon burning. However, a large fraction of the star becomes so hot that photodisintegration 
instability occurs. This uses all the energy released by previous nuclear burning stages and, 
instead of reversing the implosion, it accelerates the collapse and a proto-black hole forms 
quickly accreting mass from the dense inner part of the core. Centrifugal forces are able to 
slow down material at the equator. When the black hole mass is ~ I4OM0, a stable accretion 
disk forms with mass ~ 30Mq and accretion rates of order 1 — 10 Mq s~^. 

It is widely accepted, although not confirmed, that accretion systems of this kind rep- 
resent one of the leading models to power the relativistic fireball in GRBs (Woosley 1993). 



■^The relative error associated to this approximation is w where N is the number of individual 

bursts within the observation time 
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In order to produce a jet, the energy generated must be transported out by magnetohydro- 
dynamical (MHD) processes. This mechanism appears to be the most appropriate to very 
massive coUapsars, as other means of energy transport become inefficient as the black hole 
mass increases (Popham et al. 1999). Any quantitative prediction about the energy and col- 
limation of the magnetic-driven jet is necessarily uncertain because of the complex physics 
involved. The energy released by the jet can be written as edisk^MHoMiiskC^ where Cdisk ~ 0.11 
is the efficiency of the accretion disk for typical values of the angular momentum, eMHD ~ 0.1 
is the efficiency of the magnetic-driven explosion and Mdisk is the mass of the accretion disk 
(Fryer, Woosley & Heger 2001). This yields 6 x 10^^~^^ erg depending on the assumed 
angular momentum distribution. If beamed into 1% of the sky, this jet corresponds to an 
inferred isotropic energy E^f-o ~ 10^^ erg, 1-2 orders of magnitude larger than for ordinary 
collapsars. The jet lifetime is limited by disk-fed accretion and can be crudely estimated to 
be tj ^ 10(1 + 2;) s (Fryer, Woosley & Heger 2001). Therefore, Poplll GRBs occurring at 
redshifts ^ 10 are expected to last longer than their more recent ( ^ 5) counterparts. 

Recently, Meszaros & Waxman (2001) have proposed a mechanism for the emission of 
TeV neutrinos in collapsars. In the following, we give a brief description of the proposed 
model in order to introduce the relevant quantities for the present analysis. The interested 
reader is referred to the original paper for details. 

a) Generalities of the model 

In a generic coUapsar model, the jet propagates inside the He-core and the H-envelope 
before emerging. As the jet propagates out into the surrounding circumstellar matter, shocks 
are formed which ultimately produce 7-rays. In their model, Meszaros & Waxman (2001) 
show that as the jet is propagating through the dense He-core, it decelerates to subrelativistic 
velocities. However, beyond the He-core edge at r 10^^ cm, the density in the H-envelope 
drops down to ~ 10~^gcm~^ and the jet can accelerate to relativistic velocity. At the edge 
of the He-core the jet is capped by a termination shock and a reverse shock moves back into 
the jet, accelerating the electrons of the shocked plasma. The electrons are expected to lose 
all their energy on a very short timescale by synchrotron and inverse-Compton emission; 
the resulting radiation quickly thermalizes due to the high optical depth to an approximate 
black body spectrum whose peak energy is in the X-ray domain. Inhomogeneities in the 
jet cause internal shocks at a typical distance of ^ 10^^ cm, which is smaller than both 
the termination and the reverse shock radii (Meszaros & Waxman 2001). Internal shocks 
can accelerate protons which interact with the X-ray photons leading to pion production 
and consequently to neutrino emission. We describe this process in detail in the following 
subsections. 
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b) Target radiation 

Electrons accelerated in the reverse shock, at r 10^^ cm, lose all their energy to 
radiation which quickly thermalizes due to the high optical depth. Therefore the spectrum 
of the target photons is assumed to have a blackbody shape peaking at energy 

.,^4ry"(-^)"\eV (6) 



where E^^ — £^iso/10^^ erg is the isotropic energy of the explosion, tj^i — tj/lO s is the jet 
lifetime, = cm is the shock radius, is the Lorentz factor of the shocked jet 

plasma; all quantities are evaluated in the source rest frame. 

c) Photo-meson interaction 

Internal shocks can accelerate protons to a power law distribution which approximates 
dup/dep oc e~^, where Up is the number of protons with energy (Waxman 2001). The 
maximum proton energy is determined by equating the acceleration time, estimated as the 
Larmor radius divided by c, to the minimum between the dynamical time and the synchrotron 
cooling time. In the source rest frame this is given by, 

p5/2 ^^1/2 ^1/4 

ep,.a. ^ 9.5 X 10" ,y .t;, GeV, (7) 

where Tj = 10^'^r2.5 is the typical Lorentz factor of the jet, 6t = 10^^5t_2 s is the variability 
timescale of the injected relativistic outflow and = 10~^^b,2 is the magnetic fleld equipar- 
tition fraction. As they approach the reverse shock, high energy protons may produce both 
charged and neutral pions, with roughly equal probabilities, through the interaction with 
photons of energy e^, 

p + 7— >7r°+p p + 7 — > TT"*" + n. (8) 

In order for this reaction to take place, the proton energy must exceed the threshold energy 
for A-resonance with thermal photons, i.e. 

ep>ep,A = 0.3Tle-' GeV. (9) 

These protons will lose all their energy to pions through multiple collisions at a typical 
distance from the termination shock where the photon density is such that the photo-meson 
optical depth Tp^^ ^ 1. Charged pions will then decay into muons finally leading to neutrino 
production: 

7r+ jji^ ^ e+ + i/e + ly^. (10) 
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In a single collision, a pion is created with an average energy which is 20% of the proton 
energy. This energy is roughly evenly distributed between the final tt"^ decay products, 
yielding a neutrino energy Ci, — O.OSep. 

d) Neutrino energy spectrum 
The neutrino energy spectrum emitted by each source can be expressed as 

7.^4^. (11) 

where the neutrino energy ranges between a minimum value, ey^min = 0.05ep,A, and a maxi- 
mum value ej^^max = O.OSep^max and n^, is the number of neutrinos with energy e^,. 

Inverse Compton losses might be relevant for high energy pions and muons. Prom the 
results given in Meszaros & Waxman (2001), we derive that while pions are expected to 
decay prior to significant energy loss, muons, whose lifetime is ~ 100 times longer, may lose 
significant fraction of their energy before decaying. To be conservative we consider only those 
neutrinos which are directly produced by 7r+ decay; these carry ~ 1/4 of the pion energy. 
As a consequence, the total energy emitted in neutrinos is given by 

H 

-C/i/.tot - -^^j— 7 1— T-C/iso, (J-^j 

o mi^ep^max/ £p,minj 

where ep^min = Th^^pC^, ~ 0.4 is the fraction of the total energy in accelerated protons, 
assuming equipartition. 

Synchrotron losses are unimportant for protons within the present energy range but 
might become relevant for pions with energy greater than a threshold value above which 
pions lifetime exceeds the characteristic time for energy loss due to synchrotron emission 
(Waxman & Bahcall 1997, 1999; Rachen & Meszaros 1998). The corresponding neutrino 
energy is given by. 



e.,. = 10' 



T/i ^2.5 ^^-2 ^ 



GeV. (13) 



Since Tp^^ 3> 1, at all energies ej, < ej. ^, the v energy spectrum is simply proportional to 
that of the protons, i.e. dui^/dei, oc and we can write, 

J _ I A/" Ey^tot if ^u,m\n < < e^^s 

" 1 NE, 

Ev,tot is the total energy emitted in neutrinos with ej^^min < Ci/ < Ci/.max and H"^ — 
ln(e,/,s/e,/,inin) + 1/2[1 — (ei/,s/ei/,max)^] is a normalization factor. At neutrino energies > e^^si 
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the probability that a pion would decay before losing its energy decreases as ~ [^u/^u,s] ^ 
and the spectrum is no longer flat. 

The average energy flux, F^{E) (units GeV cm~^), observed at a given energy E and 
emitted by a source at redshift z can be expressed as, 

^ ^ dn dJ^^^[E{l + z)] _ UE{l + z)] 

^^^^^-J4^ -(^ + ") 4nr-{z) ^^^^ 

for e,y,min < E{1 + z) < e^^^max, and surface element = r'^{z)dfl, where r{z) is the comoving 
distance to the source, E{1 + z) is the emission energy, J^[E{1 + z)]/{l + z) is the redshifted 
energy spectrum and beaming effects are properly taken into account. 



4. Neutrino Diffuse Flux from PopIII GRBs 

In this Section we estimate the diffuse flux of neutrinos emitted by PopIII GRBs through 
the mechanism described above. The integrated signal observed today at energy E is 

E^^ [GeVcm-^s-^r-^] = ^ [ dR{z)F,{E), (15) 

dE 4:71 Jzf 

where dR{z) is the differential rate of all PopIII GRBs, given by eq. ^ Qj = AQj/An = 0.01 
is the fraction of sources which are beamed towards the Earth and Fy{E) is the average 



flux emitted by a single source at energy £"(1 + z) (see eq. 14). The lower limit Zf is the 
transition redshift from a top-heavy to a normal IMF marking the end of the epoch of 
VMBH formation; the value of 2;/ depends on /^^ (see eq. |l]). In computing the diffuse flux 
we are making the implicit assumption that all sources contribute the same neutrino flux to 
the overall signal, i.e. all sources have the same progenitor (SOOM©) and release the same 
amount of energy Ei^^- 

We show in Fig. || the diffuse flux obtained assuming two possible values for the equiva- 
lent isotropic energy i^iso = lO^'', 10^^ erg and Z™'^ that corresponds to the maximum source 
emission rates plotted in Fig. The experimental upper limits for present (AMANDA-BIO, 
hereafter always indicated as AMANDA-I) and forthcoming (AMANDA-II, IceCube) neu- 
trino telescopes are also shown, together with the diffuse flux from atmospheric neutrinos 
and a number of other astrophysical sources. The effect of a smaller value for ii^iso (smaller 
explosion energy or larger beaming angle) is to decrease the overall amplitude of the dif- 
fuse emission flux but at the same time to extend the range of observed energies where the 
spectrum would appear as flat. 
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As can be inferred from Fig. ^ because of the larger explosion energy, the integrated 
signal from PopIII GRBs is higher than that generated from "ordinary" GRBs, at least for a 
significant range of values for /^^ (or, alternatively Zj). The maximum flux predicted by our 
model is above the upper limit recently established by the AMANDA-I detector (Andres et 
al. 2001). This finding, together with the anticipated sensitivities of AMANDA-II (Barwick 
2001) and IceCube (Spiering 2001) shows how neutrino telescopes might enable to derive 
important constraints on the characteristic energies released by PopIII GRBs as well as on 
the IMF of the first stars. The relevant quantities that determine the amplitude (and, to a 
less significant extent, the energy range) of the diffuse flux are represented by the assumed 
explosion energy and by the PopIII GRBs IMF. As discussed in the previous Section, Eiso 
depends on the rather uncertain physics of MHD processes through the quantities eMHD, Cdisk 
and Mdisk as well as on the assumed degree of beaming. Therefore, it is equally important to 
explore a wide range of values for the explosion energy, such as 10^^~^^ erg. The aim of the 
analysis presented in the next Section is to investigate to what extent present and forthcoming 
neutrino telescopes might be able to test our scenario. In particular, if the current scenario 
is correct, important considerations regarding i^iso and /^^ might be derived even from a 
non-detection, i.e. if the signal is below the expected upper limits. This will enable to reject 
a range of models which so far appear to be untestable by any other observational facility. 



5. Experimental Predictions 

The diffuse emission from unresolved PopIII GRBs (as well as from other astrophysical 
sources) should have an energy spectrum that is much harder than the atmospheric neutrino 
spectrum (see Fig. 0). Therefore, a search for a diffuse flux can be undertaken by looking 
for an excess of events at large energies (see e.g. Andres et al. 2001). Unfortunately, the 
energy resolution of a neutrino telescope is rather poor. Indeed, the energy of the event is 
roughly estimated counting the number of optical modules triggered by Cerenkov radiation 
from muons that are produced in neutrino-nucleon interactions in the ice surrounding the 
detector or in the bedrock below (Andres et al. 2000). Thus, in the following, we consider 
two possible energy intervals: [lO'' — 10^] GeV and [10^ — 10^] GeV that we indicate as 
low- (L) and high- (H) energy bands, respectively. Within these energy bands, the diffuse 
emission from PopIII GRBs is higher than the predicted experimental sensitivities, at least 
for a range of model parameters (see Fig. ^). The atmospheric neutrino emission in the 
L-band is the primary source of contamination as the contribution from other astrophysical 
sources is well below the experimental sensitivity (see Fig. |^). 

In Fig. I we show the differential number of neutrinos, dN/dE, predicted by our model 
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assuming the same model parameters as in Fig. ^ i.e. Eiso = 10^^ erg and Z™"^. The same 
quantity corresponding to the atmospheric neutrino diffuse ffux and to the predicted upper 
hmits for AMANDA-I (upper panel), AMANDA-II (medium panel) and IceCube (lower 
panel) are shown both in the L- and the H-bands. The shaded areas indicate the range 
of models corresponding to the same E[so but to different values for /^^ that might lead 
to signal detection or, equivalently, that might be rejected by a non-detection. As can be 
inferred from the Figure, the upper limit set by AMANDA-I (Andres et al. 2001) can be 
used to reject a region of model parameters: if i^iso = 10^^ erg, than /^^ must be > 0.27 and 
the transition to a normal IMF must have occurred at redshifts Zf > 9.2. 

The building of km-scale detectors will significantly improve the sensitivities in both 
bands, even though the L-band appears to be largely dominated by the atmospheric neutrino 
diffuse flux. Thus, AMANDA-II and IceCube are expected to show equivalent performances 
in this channel. Conversely, because of the softer energy spectrum, the atmospheric neutrino 
signal drops in the H-band and IceCube is expected to significantly improve the experimen- 
tal constraints that might be achieved by AMANDA-II, as can be seen by comparing the 
extension of the shaded regions in Fig. ^. 

The effect of variations in the values of Eiso and f^-y can also be explored. To make a 
systematic analysis of the parameter space, for each model [identified by a point in the (-Eiso, 
/^^) plane] we compute the quantity -Rl (-Rh) defined as the total number of neutrinos in the 
L-band (H-band) normalized to the largest between the number of atmospheric neutrinos 
and the experimental upper limits for AMANDA-I, AMANDA-II and IceCube. 

Figs. ^ and | show maps of i?L (upper panels) and i?H (lower panels) for AMANDA-II 
and IceCube. The solid lines mark the limits of the regions (right to the curves) where these 
detectors would be able to detect a count excess. The dashed lines show the current limit set 
by AMANDA-I, which corresponds to explosion energies Eiso in the range [6.4x lO^''— 10^^] erg 
and to f^^ < [8 X 10"^ — 0.27] depending on Eiso- 

Thus, according to the present observational constraints, if PopIII GRBs had occurred 
with energies in the above range, than the fraction of SN^^to VMBHs must have been higher 
than the above limit {i.e. the transition from a top-heavy to a normal IMF must have 
occurred at redshifts zj > [5.4 — 9.2] depending on Eiso) in order for the overall signal to be 
below the experimental upper limit. 

In the L-band, AMANDA-II and IceCube would be able to extend the region excluded 
by AMANDA-I down to energies 3.6 x 10^^ erg and up to f-^y < 0.43 (transition redshifts 
Zf ^ 9.75). Moreover, the enhanced sensitivities of km-scale detectors open the possibility to 
test our scenario in almost all range of plausible explosion energies {Eiso > 7.6 x 10^^ erg) and 
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to derive fundamental insights on the IMF of the first stars, with /^^ ranging from 8 x 10 ^ 
up to 0.98. 

Thus, unless MHD efficiencies are so low to significantly reduce the energy released by 
PopIII GRBs, the anticipated sensitivities of AMANDA-II and IceCube will very hkely lead 
to a direct detection of PopIII GRBs in at least one of the two energy bands. The above 
range of values on corresponds to two opposite scenarios for the evolution of the first 
stars: if f^^ ^ 8 x 10^^ than approximately 1 SN^^ forms every 120 VMBHs (hence PopIII 
GRBs). This means that a very small amount of mass has been processed into metals and 
that star formation proceeds according to a top-heavy IMF down to redshifts Zf ^ 5.4. Since 
the bulk of the mass went in PopIII progenitors with masses > 26OM0, we expect a large 
number of VMBHs remnants to be present today. The resulting critical density in VMBHs 
has been estimated to be high enough to account for the entire baryonic dark matter halos 
of present-day galaxies (see Schneider et al. 2001 for a thorough discussion of these issues). 
Conversely, if f^-y 0.98, than approximately 1 VMBH forms every 75 SN^^. In this case, 
a large fraction of mass in PopIII stars has been processed into metals and transition to a 
normal IMF occurs at much higher redshift, at Zf ^ 10. However, the critical density in 
VMBHs formed at these early stages is still comparable to that of super massive black holes 
observed in the nuclei of present-day galaxies (Schneider et al. 2001). 



6. Implications for PopIII Gamma Ray Bursts 

At present it is not clear if the relativistic jets can break-out of the stellar envelope. In 
PopIII GRBs progenitors, jet driving depends mainly upon MHD processes. For this reason, 
such jets tend to be relatively cold, in the sense that their thermal energy is small compared 
to their kinetic energy (MacFadyen, Woosley & Heger 2001) and hence they are efficiently 
pressure-coUimated by the He-core through which they propagate. However, the internal 
structure of PopIII stars is different from that of ordinary coUapsars, being characterized 
by larger He-core densities and more extended H-envelopes. Under these conditions, it is 
hkely that jet staUing might occur before break-out if the crossing time is longer than the 
jet lifetime. Since the jet becomes relativistic near the outer radius of the He-core, a more 
extended H-envelope does not affect its crossing time, which is substantially determined by 
the fly time through the He-core (Mcszaros & Rees 2001). However, the effect of rotation in 
PopIII GRBs is to create a funnel along the rotation axis within which the density becomes 
substantially lower than the equatorial one. Therefore, at least in some cases, this might 
lead to a break-out. If this is the case, GRBs are a necessary endproduct of such process; if 
instead the jet is chocked the neutrino signal will nevertheless be emitted. In the previous 
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Sections we have already analyzed in detail the latter scenario. Here we briefly analyze the 
imphcations descending from a successful jet break-out, leading to a population of very high 
redshift GRBs. What are the distinctive signatures that would enable us to pin-point PopIII 
GRBs among observed ones? 

Probably the most obvious such features are the following: i) depending on the proper- 
ties of the accretion disk and on the uncertain interaction of the jet with surrounding matter, 
the PopIII GRBs should probably on average last longer, due to cosmological time dilation; 
a) the peak of emission, that in the rest frame is in 7-rays, would be shifted into X-rays; in) 
their optical afterglow should be heavily absorbed by the intervening intergalactic medium 
(particularly if they occur before reionization) ; iv) iron lines should be very weak or totally 
absent, because nucleosynthetic products of progenitor stars are swallowed by the VMBH 
and because the surrounding gas is presumably of primordial composition. Have events with 
similar characteristics already been observed? 

Intriguingly, BeppoSAX has revealed the existence of a new class of events, the so-called 
X-ray flashes or Fast X-ray Transient (FXTs), which show the bulk of emission in X-rays 
(Heise et al. 2001). They are bright X-ray sources, with peak energies in the 2 — 10 keV 
band and durations between 10 — 200 s, which are not triggered and not detected in the 
7-ray range 40 — 700 keV. A subsample of 10 FXTs have been observed in the lowest energy 
channels of the BATSE instrument (Kippen et al. 2001). This has enabled to directly 
compare the gamma-ray properties of FXTs to those of "ordinary" GRBs. The preliminary 
result is that FXTs appear to be consistent with an extrapolation of GRB behaviour, showing 
softer spectra (based on their fluence hardness relation) and duration all in excess of 10 s. 
While the limited sample does not allow to draw statistically solid conclusions, none of 
the observational evidences is in apparent conflict with their nature being PopIII GRBs. 
Moreover, peak energies in the observed range (2 — 10 keV) are naturally predicted in PopIII 
GRBs not only as a consequence of higher emission redshifts. Indeed, gamma-ray emission 
will occur at larger radius because of the more extended H-envelope, leading to a lower value 
for the emitted peak energy (^ 200 keV, Guetta, Spada & Waxman 2001). In addition, 
requirement in) above is very likely fulfilled. In fact, a significant number (about 40% 
of GRBs for which fast follow-up observations were carried out) of GRBs do not show an 
optical counterpart; for this reason they are dubbed as GHOST (GRB Hiding Optical Source 
Transient). It is quite interesting to note that all FXTs fall into the GHOST category^ 



''As we are writing, an X-ray rich GRB (GRB011211) has been detected by BeppoSax WFCl. FoUow- 
up observations have identified the optical afterglow and have measured a redshift z — 2.1 (Piro, private 
comunication) . However, a detailed analysis of the BeppoSAX GRBM data shows a very shallow long event 
in the 40-700 keV band with total duration of about 270 s, similar to that measured in X-rays (2-28 keV) by 
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The above scenario cannot yet be considered more than a promising possibihty and 
there are alternative explanations. For example, the extension of GRBs to the low peak 
energy regime may be a natural consequence of the fireball model. As pointed out by 
Guetta, Spada & Waxman (2001), FXTs could be also produced by relativistic winds where 
the minimum Lorentz factor is smaller than that required to produce "normal" GRBs. The 
region of parameter space leading to FXTs in their model is rather small; in addition, they do 
not make any specific prediction on the burst duration. In this framework, the failed optical 
detection can be caused by dust extinction within the host system, although this explanation 
is still debated (Lazzati, Covino & Ghisellini 2001; Djorgovski et al. 2001, Ghisellini 2001). 
Prompt infrared observations would be able to tackle this issue. 

The number of undetected FXTs is estimated to be very large. The untriggered BATSE 
catalogue suggests an all-sky rate of ~ 400 yr^^ (Stern et al. 2000). If we maintain the 
above speculation that FXTs are indeed successful {i.e. non chocked) PopIII GRBs, and if 
a neutrino precursor from a single event could be measured in coincidence with an X-ray 
flash, this would provide a direct probe of the nature of X-ray flashes and on the relative 
number of chocked and successful GRBs. An estimate of the likelihood of such an event 
can be given as follows: the average number of upward muon detections per unit area in a 
neutrino telescope during an individual burst is 



where P^^ = 1.3 x IQ-^E/TeVy {P = 2 if E < ITeV and (3 = 1 if E > 1 TeV) is 



the probability that a muon neutrino will produce an upward moving muon in a detector 
(Gaisser, Halzen & Stanev 1995) and E^{E) is the average neutrino flux emitted by a single 



as IceCube) and a source energy, the number of muon detections per burst is A^"^ = 0.03 — 
0.003 for emission redshift 5.4 < z < 30. This signal is signiflcantly higher than both the 
atmospheric neutrino background and the cumulative PopIII GRBs signal at comparable 
energies. Since the typical angular resolution of planned neutrino telescopes is 6 ^ 1°, 
a diffuse flux will produce a negligible average number of muons around the direction of 
arrival from a point source within the duration of the burst. Assuming the above mentioned 



WFCl. Given the redshift, the (isotropic) 7-ray energy is 6.3 x 10^^ erg (Frontera et al. GCN 1215). This 
prehminary analysis shows that the X/7 fluence ratio is similar to that found in a number of ordinary GRBs 
detected by BeppoSAX. The distinguishing features of GRB011211 are its long duration (the longest event 
localized with BeppoSAX) and its faintness both in X- and 7-rays. Therefore, the analysis leads to conclude 
that this event cannot be classified as a FXT. 
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rate of 400 yr^^ for FXTs, one could expect a few neutrino bursts/yr with A*"^ > 1. This 
result is obtained assuming that all FXTs are successful PopIII GRBs beamed towards the 
Earth, which for a beaming factor of 1% and values for /^^ in the detectable region, implies 
a ratio of successful to chocked GRBs of less than 5%. 



7. Summary of the Results 

The main results of this paper can be summarized as follows: 

• The total rate of PopIII GRBs produced by massive, first stars can be as high as 
4 X 10^ yr^^ for the most favourable parameter choice and if the jets are beamed into 1% of 
the sky; for some model parameters, the event duty cycle is found to be smaller than unity 
implying that the diffuse neutrino signal might be characterized by a sequence of individual 
bursts. 

• High energy neutrinos from PopIII GRBs could dominate (above other plausible neu- 
trino sources) the flux in the L [10^ — 10^] GeV and H [10^ — 10^] GeV energy bands which 
can be explored by present (AMANDA-I) and forthcoming (AMANDA-II, IceCube) neutrino 
telescopes (Fig. 

• Constraints on /^^ [i.e. on the mass fraction of PopIII stars which ends up as SN^^) 
from these experiments will also provide a tremendous insight into the IMF of the first stars, 
and in particular, on the redshift Zf, which marks the metallicity-driven transition from a 
top-heavy IMF to a the current, standard [i.e. Salpeter-like) one. 

• For explosion energies in the range Eiso = 6.4 x lO^'' — 10^^ erg, the current upper 
limit established by AMANDA-I in the L- and H-bands can be used to reject a region of the 
model parameter space, namely /^^ < [8 x 10^^ — 0.27] corresponding to transition redshifts 
Zf < [5.4 — 9.2]. Thus, either PopIII GRBs are characterized by smaller explosion energies or 
massive primordial star formation occurred only at very high redshifts and metals released 
by SN^-y enriched the IGM to the critical level < Z >= IO'^^Zq before redshift 9.2. 

• The enhanced sensitivities of forthcoming experiments such as AMANDA-II and, in 
particular, km-scale detectors such as IceCube, might be able to discriminate between the 
two above possibilities: observations in the H-band will explore a wide range of explosion 
energies, i^iso = 7.6 x 10^^ — 10^^ erg, values for /^^ in the range 8 x 10~^ < /^^ < 0.98 and 
transition redshifts 5.4 < < 10. 

• For a subset of model parameters, neutrino emission from individual bursts might be 
detected with km-scale detectors in the L-band over roughly a year of observations. To be 
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statistically significant, single burst detection should occurr in coincidence with observations 
of the GRB itself, if the jet is not chocked and can successfully propagate out of the stellar 
envelope. 

• We have speculated that Poplll GRBs, if they are not chocked, could be associated 
with a new class of events, the Fast X-ray Transient (FXTs), which are bright X-ray sources, 
with peak energies in the 2 — 10 keV band and durations between 10 — 200 s, which are not 
triggered and not detected in the 7-ray range 40 — 700 keV. An additional hint that they 
could be high redshift sources comes from the fact that all FXTs observed so far fall into the 
GHOST (GRB Hiding Optical Source Transient) category, an occurrence highly suggestive 
on the absorbing effects of the predominantly neutral intergalactic medium present at epochs 
prior to cosmic reionization. 

We are grateful to B. Ciardi, C. di Stefano, K. Mannheim, G. Riccobene, C. Spiering 
for stimulating discussions and useful information. We particularly acknowledge L. Piro and 
E. Waxman for useful conversations and ideas which greatly improved this work. This work 
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Fig. 1. — Predicted total rate of Pop III GRBs as a function of redshift is shown (solid lines) for two 
different values of f^^ (hence Zf, see text). For the same values, we plot (dashed lines) the duty cycle of the 
signal assuming the duration of each burst to be 10(1 + z) s. 
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Fig. 2. — The diffuse neutrino flux from different sources is compared to the sensitivities of present and 
forthcoming telescopes. The diffuse emission from Pop III GRBs is computed assuming two possibile values 
for the explosion energy, E'iso = 10^'', 10*^^ erg (solid curves labeled 54 and 56 respectively) and /™" that 
corresponds to the maximum source rate (see Fig. The long-dashed lines represent the diffuse emission 
of atmospheric neutrinos in the horizontal (upper boundary) and vertical (lower boundary) direction (Lipari 
1993). The other models correspond to the diffuse emission from "ordinary" GRBs and GRBs afterglows 
(Waxman & Bahcall 1997), to a model for photo-meson interaction in blazar jets producing ultra- high energy 
cosmic rays through neutron escape and to decaying gauge bosons created at topological defects (see Learned 
& Mannheim 2000 and references therein). 
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Fig. 3. — The number of neutrinos per energy interval as a function of the observed energy predicted by 
our reference model (-Eiso — 10^^ erg) assuming (maximum source rate). The signal is compared to the 
atmospheric neutrino diffuse emission (dashed line, corresponding to the upper bound in Fig. |^) and to the 
sensitivities of AMANDA-I (upper panel), AMANDA-II (medium panel) and IceCube (lower panel). The 
shaded areas identify the regions of parameter space (corresponding to the same i^iso but to different f^^) 
that would appear as an excess of events in the low and high-energy bands. 
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Fig. 4. — Upper panel: values of i?L i-e. the ratio of predicted number of neutrinos in the low energy band 
(L: 10'' - 10^ GeV) normahzed to the predicted AMANDA-II sensitivity in the same band, in the (-Eiso, f-yy) 
plane. The color bar shows i?L in Log scale. The solid line identifies the region (right to the curve) where 
i?L > 1; the region rightmost of the dashed line is excluded by the AMANDA-I current limit. Lower panel: 
same as above but for i?H, i-e. the ratio of predicted number of neutrinos in the high energy band (H: 
10^ - 10'^ Gev) to AMANDA-II sensitivity in the same band. 
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Fig. 5. — Same as Fig. | but for IceCube. 



